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ABSTRACT 

We report on the discovery of ETHOS 1 (PN G068. 1+11.0), the first spectroscopically 
confirmed planetary nebula (PN) from a survey of the SuperCOSMOS Science Archive 
for high-latitude PNe. ETHOS 1 stands out as one of the few PNe to have both 
polar outflows (jets) travelling at 120 ± 10 km/s and a close binary central star. The 
lightcurve observed with the Mercator telescope reveals an orbital period of 0.535 days 
and an extremely large amplitude (0.816 mag) due to irradiation of the companion 
by a very hot pre-white dwarf. ETHOS 1 further strengthens the long suspected link 
between binary central stars of planetary nebulae (CSPN) and jets. INT IDS and VLT 
FORS spectroscopy of the CSPN reveals weak N III, C III and C IV emission lines 
seen in other close binary CSPN and suggests many CSPN with these weak emission 
lines are misclassified close binaries. We present VLT FORS imaging and Manchester 
Echelle Spectrometer long slit observations from which a kinematic model of the nebula 
is built. An unusual combination of bipolar outflows and a spherical nebula conspire to 
produce an A-shaped appearance. The kinematic age of the jets (1750 ± 250 yrs/kpc) 
are found to be older than the inner nebula (900±100 yrs/kpc) consistent with previous 
studies of similar PNe. Emission line ratios of the jets are found to be consistent 
with reverse-shock models for fast low-ionisation emitting regions (FLIERS) in PNe. 
Further large-scale surveys for close binary CSPN will be required to securely establish 
whether FLIERS are launched by close binaries. 

Key words: planetary nebulae: individual: PN G068. 1+11.0 - binaries: close - stars: 
winds, outflows - ISM: jets and outflows 



1 INTRODUCTION 

In the protracted debate on the shaping mechanisms of PNe 
(Balick & Frank 2002) appropriate solutions are sought to 
explain both the dominant nebula morphology (e.g. spher- 
ical, elliptical, bipolar; Balick 1987) and the accompanying 
collimated outflows or 'jets'. Single star evolution readily ex- 



* Based on observations made with the Flemish Mercator and 
Isaac Newton Telescopes of the Observatorio del Roque de Los 
Muchachos and the Very Large Telescope at Paranal Observatory 
under programs 083.D-0654(A) and 085.D-0629(A). 
f E-mail: b.miszalski@herts.ac.uk 



plains spherical PNe and potentially elliptical PNe if there 
is interaction with the interstellar medium (ISM; Villaver 
et al. 2003; Wareing et al. 2007). Highly axisymmetric or 
bipolar nebulae (Corradi & Schwarz 1995) have been mod- 
elled with some success including the generalised interact- 
ing stellar wind model (GISW; Kwok, Purton & Fitzgerald 
1978; Kwok 2000) and a combination of stellar rotation and 
magnetic fields (e.g. Garcia-Segura 1997). There are how- 
ever fundamental limitations in both of these models. The 
GISW model depends upon an assumed density contrast 
to produce bipolar PNe (e.g. a dusty torus) and the latter 
cannot operate without an additional supply of angular mo- 
mentum (e.g. a binary companion; Soker 2006; Nordhaus, 
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Blackman & Frank 2007). While there are multiple advan- 
tages to a binary explanation over existing theories, we are 
only beginning to understand the prevalence and impact of 
binarity in PNe (De Marco 2009; Miszalski et al. 2009a, 
2009b, 2010). 

Decisive observational evidence is required to advance 
the shaping debate which is imbalanced towards theoretical 
models. PNe with dusty disks seem to support the GISW 
model though there appears to be a strong dependence on a 
binary companion for their formation (Chesneau 2010). Sim- 
ilarly, magnetic fields are difficult to observe directly in PNe 
(Jordan et al. 2005; Sabin et al. 2007), yet they appear to be 
more influential as part of a binary-driven dynamos (Nord- 
haus & Blackman 2006; Nordhaus et al. 2007). Evidence 
for close binary central stars of PNe (CSPN) that passed 
through a common-envelope (CE) phase is now firmly es- 
tablished and gaining considerable momentum with at least 
40 now known (Miszalski et al. 2009a, 2010; De Marco et al. 
2008). These make up at least 17±5% of all CSPN (Miszalski 
et al. 2009a; Bond 2000). 

With a larger sample of 30 post-CE nebulae Miszalski 
et al. (2009b) found at least 30%, perhaps 60-70%, of post- 
CE binaries had bipolar nebulae suggesting the CE phase 
preferentially shapes bipolar nebulae. This is a substantial 
improvement over the study of Bond & Livio (1990) whose 
sample of around a dozen close binaries showed no clear 
morphological trends. High resolution kinematics of nebulae 
around post-CE CSPN have shown a tendency for matched 
orbital and nebulae inclinations that strengthen the connec- 
tion between the binary and its nebula (e.g. Mitchell et al. 
1997; Jones et al. 2010). Further progress in this aspect will 
require a statistically significant sample of post-CE nebu- 
lae to be identified and studied in detail to reveal trends 
bestowed upon the nebula during the CE phase (Miszalski 
et al. 2010). Miszalski et al. (2009b) made the first steps 
in this direction by finding a propensity for low-ionisation 
structures (Gongalves et al. 2001) amongst post-CE nebulae 
either in the orbital plane as a ring (e.g. Sab 41, Miszalski 
et al. 2009b; the Necklace, Corradi et al. 2010) and in the 
polar direction as jets (e.g. A 63, Mitchell et al. 2007). 

Jets occur in a wide variety of astrophysical objects in- 
cluding PNe in which they are the least understood (Livio 
1997). The connection between jets and binarity has long 
been suspected (Soker & Livio 1994) but never proven given 
the paucity of known PNe with jets and binary CSPN. Jets 
in PNe typically have low outflow velocities (~100 km/s), 
but higher velocities of 200-300 km/s are not uncommon, 
e.g. M 1-16 (Schwarz 1992; Corradi & Schwarz 1993), NGC 
6337 (Corradi et al. 2000) and NGC 2392 (Gieseking, Becker 
& Solf 1985), and can even reach 630 km/s in the case of 
MyCn 18 (O'Connor et al. 2000). There is no strong evi- 
dence for highly collimated outflows (e.g. Herbig-Haro-like 
outflows) in genuine PNe and Frew & Parker (2010) sug- 
gest the few objects with the most collimated outflows are 
instead related to B[e] stars. The most striking example is 
Hen 2-90 (Sahai & Nyman 2000) which was later reclassified 
as a B[e] star (Kraus et al. 2005; Frew & Parker 2010). 

Instead, jets in PNe manifest as point-symmetric 
'corkscrew'-like outflows (Lopez et al. 1993; Vazquez et al. 
2008) or as pairs of opposing knots often called fast low- 
ionisation emitting regions (FLIERS; see e.g. Balick et al. 
1987, 1993; Dopita 1997) with notable examples being Hb 4 



(Lopez et al. 1997; Harrington & Borkowski 2000; Miszalski 
et al. 2009b), NGC 6337 (Corradi et al. 2000; Garcfa-Di'az 
et al. 2009), A 63 (Mitchell et al. 2007), IC 4634 (Guerrero 
et al. 2008), IC 4673 (Kovacevic et al. 2010), Sab 41 (Mis- 
zalski et al. 2009b) and the Necklace (Corradi et al. 2010). 
Further examples are listed by Gongalves et al. (2001) and 
Harrington & Borkowski (2000). 

Point-symmetric jets have received the most attention 
in the literature with models consisting of an episodic pro- 
cessing jet powered by an accretion disk in a wide binary 
with P or b~a few years (e.g. Raga, Canto & Biro 1993; Cliffe 
et al. 1995; Haro-Corzo et al. 2009; Raga et al. 2009). Pre- 
cession helps smear out or widen the jet and may also ex- 
plain the often 'bent' position of jets with respect to the 
presumed major axis (e.g. Hb 4 and Sab 41). Some jets ap- 
pear in two pairs (e.g. NGC 6337) and these may have arose 
from multiple ejections (e.g. Nordhaus & Blackman 2006). 
A similar scenario involves a short-lived 'wobbling' accre- 
tion disk (Livio & Pringle 1996, 1997; Icke 2003) and mod- 
els with magnetic fields have also successfully been applied 
to the problem (Garcfa-Segura 1997). Dynamical lifetimes 
of jets suggest they are ejected well before the main nebula 
(Mitchell et al. 2007; Corradi et al. 2010), i.e. before the CE 
is ejected if a close binary is present, and this is consistent 
with the short-lived nature of the modelled accretion disks. 

The most constructive way to establish the link between 
jets and binaries is to look for binary CSPN in a large sam- 
ple of PNe with jets. This is one of the major aims be- 
hind our ongoing program to look for close binary CSPN 
via the photometric monitoring technique (Miszalski et al. 
2010). We are also targeting bipolar PNe and through non- 
detections we will also be able to discern whether larger 
orbital periods produce point-symmetric jets rather than 
FLIERS. The program has been very successful so far with 
at least five close binary CSPN found using the 1.2-m Flem- 
ish Mercator telescope. These include the Necklace and its 
Supernova 1987A-like ring (Corradi et al. 2010), the Type 
la supernova candidate double-degenerate nucleus of Hen 2- 
428 (S ant ander- Garcia et al. 2010), and the subject of this 
paper, ETHOS 1. 

This paper is structured as follows. Sec.[2]mtroduces the 
new survey responsible for the discovery of ETHOS 1 and 
the spectroscopic confirmation of its PN status. Spectro- 
scopic and photometric evidence for binarity are presented 
in Sec. [3] The nebula is examined in Sec. U with narrow-band 
images, high-resolution spectroscopy. Sec. 14.21 examines di- 
agnostic emission line ratios of the jets and we conclude in 
Sec. [5] 



2 DISCOVERY AND PN CONFIRMATION 

2.1 The Extremely Turquoise Halo Object Survey 
(ETHOS) 

Beyond the confines of Galactic Plane Ha surveys (Parker 
et al. 2005; Drew et al. 2005) there are very few PNe cur- 
rently known as less than 9% of ~2800 Galactic PNe are 
located at latitudes \b\ > 10°. Most of these were found 
from visual inspection of photographic broadband sky sur- 
vey plates (e.g. Kohoutek 1963; Abell 1966; Longmore 1977; 
Weinberger 1977; Longmore & Tritton 1980; Hartl & Trit- 
ton 1985). With the digitisation of the same material it is 
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now possible to extend this work to lower surface bright- 
nesses and to systematically search larger areas more effi- 
ciently (Jacoby et al. 2010). The MASH-II project (Miszal- 
ski et al. 2008) demonstrated that visual searches for PNe 
are largely insensitive to compact or star-like PNe. Several 
compact PNe were found from catalogue queries of the Su- 
perCOSMOS Ha Survey (Parker et al. 2005) including the 
MASH-II PN MPA 1327-6031 that lies just near the outer 
lobe of the bipolar MASH PN PHR 1327-6032 (Parker et 
al. 2006). The number of distant (compact) Halo PNe is 
therefore essentially unknown. Further discoveries are es- 
sential to determine whether Halo PNe can only be formed 
via close binary evolution (Bond 2000) as there is currently 
only a small amount of supporting evidence (e.g. Pena & 
Ruiz 1998; Tovmassian et al. 2010) 

The Extremely Turquoise Halo Object Survey (ETHOS, 
Miszalski et al. in prep) is a new survey designed to find 
high-latitude compact PNe that would ordinarily be missed 
by visual searches. ETHOS is essentially an extension to 
broadband colours of the MASH-II catalogue search ap- 
plied to the whole sky (excluding the Ha surveyed regions) . 
Currently four candidate PNe have been spectroscopically 
confirmed out of a few dozen candidates identified solely 
from queries of the SuperCOSMOS Science Archive (SSA, 
Hambly et al. 2004), a convenient SQL interface to the en- 
tire database of catalogue data generated during the digi- 
tisation of SuperCOSMOS Sky Survey data (Hambly et al. 
2001b). Candidates are then visualised as for MASH-II, but 
this time replacing the primary colour-composite Ha (red), 
SR (green) and Bj (blue) image with an In (red), Rf 
(green) and Bj (blue) image (Reid et al. 1991; Hambly et 
al. 2001a). The IRB image is particularly effective in dis- 
tinguishing PNe apart from other objects since they often 
appear with a strong turquoise hue. Figure [T] best shows the 
turquoise false-colour in the discovery image of ETHOS 1, 
the first PN candidate identified during the design phase 
of ETHOS (Miszalski 2009). The object was selected based 
on the colours of the inner nebula but the presence of faint 
extensions either side made it stand out amongst the other 
candidates. Further details of the survey and initial results 
will be presented elsewhere (Miszalski et al. in prep). 

2.2 Spectroscopic confirmation of ETHOS 1 

ETHOS 1 was observed during Isaac Newton Telescope 
(INT) service time with the Intermediate Dispersion Spec- 
trograph (IDS) on 10 March 2009. The R300V grating was 
used at a central wavelength of 5400 A to provide wave- 
length coverage from 3300-8000 A at a dispersion of 1.87 
A/pixel. The spectrograph slit was placed along the jets at 
a position angle (PA) of 149° and its width of 1.0" gave a res- 
olution of 5 A (FWHM) at 5400 A. One 30 minute exposure 
was taken at an airmass of 1.51 followed by a 30 s exposure 
of the spectrophotometric standard star HD 192281 at an 
airmass of 1.72 with an 8" slit. Standard reductions were 
performed with IRAF and four different extractions corre- 
sponding to different nebula regions were made using the 
IR AF task APALL. The stellar continuum of the CSPN was 
traced to guide the extractions. Aperture sizes used were 6 
pixels for the CSPN (2.4"), 54 pixels for the inner nebula 
(21.6", which includes the CSPN), 19 pixels for the SE jet 
(7.6") and 23 pixels for the NW jet (9.2"). The jet extrac- 



Figure 1. Discovery image of ETHOS 1 made from ijv (red), 
R F (green) and Bj (blue) POSS-II SuperCOSMOS Sky Survey 
images (Reid et al. 1991; Hambly et al. 2001a). The turquoise 
hue is typical of PNe in this particular colour-composite image 
combination. North is to top and East to left in this lxl arcmin 2 
image. 



Table 1. Basic properties of ETHOS 1. 



PN G 068.1+11.0 

RA (J2000) 19 h 16 m 31.5 s 

Declination (J2000) +36°09'47.9" 

Diameter (main nebula) 19.4" 

Jet length (tip-to-tip) 62.6" 

c(H3) 0.18 

Heliocentric Radial Velocity —20 km/s 



tion widths were measured from the jet tips inward towards 
the CSPN. The extractions were then flux calibrated in the 
standard fashion with the aforementioned standard. 

Figure[2]shows the IDS spectra that confirm the PN na- 
ture of ETHOS 1 (PN G068.1+H.0) whose basic properties 
are given in Tab. [T] We measured c, the logarithmic extinc- 
tion at H/3, to be 0.18 using the Howarth (1983) reddening 
law and T e = 17700 K (Sec. IP]) . The high ratio of He II 
A4686/H/3 > 1.1 in the main nebula is typical of very high 
excitation PNe at high Galactic latitudes in which all the he- 
lium is doubly ionised giving a mass-bounded, optically thin 
nebula (Kaler 1981; Kaler & Jacoby 1989). Similar emission 
line intensities are also found in the post-CE PNe K 1-2 
(Exter et al. 2003) and the Necklace (Corradi et al. 2010), 
but in our case we do not find any low-ionisation species in 
the main nebula perhaps suggesting a slightly higher CSPN 
temperature. The only low-ionisation species are found in 
the tips of the jets which are caused by shocks (Sec. 14. 2[) . 
Besides the [Ne V] and [Ar V] emission found close to the 
hot CSPN, the CSPN spectrum also shows C III A4650, C 
IV A5801, A5812 and A7726 as the signature of an irradiated 
close companion (Sec. 13. 1) . 
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Figure 2. INT IDS spectroscopy of ETHOS 1. No nebular subtraction was made for the CSPN spectrum. 



3 THE BINARY CENTRAL STAR 

3.1 An irradiated spectroscopic binary 

ETHOS 1 was also observed with the FORS2 instrument 
(Appenzeller et al. 1998) under the VLT visitor mode pro- 
gram 083.D-0654(A) on 14 July 2009. The objective of this 
program was to measure the radial velocity curves of binaries 
discovered by Miszalski et al. (2009a), however inclement 
weather forced the execution of a backup program in which 
ETHOS 1 was included. The blue optimised E2V detector 
from FORS1 was used in combination with the 1200g grism 
to give continuous wavelength coverage from 4058-5556 A. 
The CCD was readout with 2x2 binning to give a disper- 
sion of 0.72 A/pixel, a spatial resolution of 0.252"/pixel and 
a resolution measured from arc lines of 1.57 A (FWHM) at 
4800 A. 

A 40 minute spectrum of ETHOS 1 was obtained at an 
airmass of 2.04 and the seeing varied between 1.0-1.3"under 
thin-thick cloud (relative flux rms up to 0.07) limiting the 
SNR in the continuum to -40 near He II A4686. The 0.7" slit 
was placed at PA = 144° to cover the jets but because of the 
gap between FORS chips not all of the NW jet was observed. 
The spectroscopic data were processed using the ESO FORS 
pipeline and the IRAF task APALL was used to make three 
extractions using the CSPN continuum as a reference. These 
were 8 pixels for the CSPN (2.0"), 84 pixels for the inner 
nebula (21.2", which includes the CSPN) and 70 pixels for 



the SE jet (17.6"). The jet extraction widths were measured 
from the jet tips inward towards the CSPN. Each spectrum 
was flux calibrated using a 300s exposure of the CSPN of 
NGC 7293 (Oke 1990) obtained at the end of the same night 
at an airmass of 1.27 and a slit width of 0.7". The final 
spectra are depicted in Fig. [3] 

Features identified in the IDS and FORS CSPN spectra 
are typical of an irradiated close binary CSPN (Pollacco & 
Bell 1993, 1994). The hot primary is revealed by He II A4686, 
A5412 in absorption and the primary is irradiating a main- 
sequence companion to produce N III (A4634.14, A4640.64), 
C III (A4647.42, A4650.25) and C IV (A4658.30, A5801.33, 
5811.98, A7726.20) emission lines. These lines originate from 
the irradiated zone on the side of the companion facing the 
primary and therefore serve as a proxy for the motion of 
the secondary. However, it should be noted that these lines 
sample the centre of light and not the centre of mass of the 
secondary (e.g. Exter et al. 2003). 

To measure any shifts present we used the [O III] and 
H/3 emission lines of our FORS spectrum of the inner nebula 
as a reference. The systemic heliocentric velocity of the inner 
nebula was set to be V aya = — 20 km/s from higher resolution 
spectroscopy (Sec.QTT]). The N III A4640.64, C III A4647.42, 
A4650.25, and C IV A4658.30 emission lines gave an average 
V^ec = —122.0 ± 4.6 km/s. Estimating Vpri is much more 
difficult given the prominent nebular He II emission but the 
He II A4686 absorption line is present (albeit at low S/N, 
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Figure 3. VLT FORS spectroscopy of ETHOS 1. No nebular subtraction was made for the rectified CSPN spectrum. 



Fig. [3]) and redshifted compared to the nebula velocity. The 
exact line profile centre is difficult to establish given the 
nebular emission and the chance that the line core may be 
partially filled in by irradiated emission. The large uncer- 
tainty means we can only estimate Vp r i ~ 105 ± 80 km/s. 
Our single epoch spectrum therefore confirms the presence of 
a spectroscopic binary with both components well separated 
from each other and the nebula. Given the large uncertainty 
in V^ r i and the lack of a well-sampled radial velocity curve 
it is premature to estimate the masses of each component 
at this stage. 

An important corollary is that the presence of the N 
III, C III and C IV emission lines alone would ordinarily be 
sufficient grounds for a weak emission line star or wels clas- 
sification (Tylenda et al. 1993; Marcolino & Araiijo 2003). 
The wels are a heterogeneous class of CSPN that are poorly 
defined (Fogel, De Marco & Jacoby 2003) and do not fit well- 
established classification schemes (Mendez 1991; Crowther 
et al. 1998; Acker & Neiner 2003). According to Tylenda 
et al. (1993), wels are identified by (i) C IV A5801, 5812 
weaker and narrower than in [WC] CSPN, (ii) very weak or 
undetectable C III A5696, and (iii) the emission complex at 
A 4650. However, this is exactly what is observed in many 
close binary CSPN with the best examples being Abell 46 
(Pollacco & Bell 1994), ETHOS 1 (Fig.EJ and the Necklace 
(Corradi et al. 2010). Many wels are therefore likely to be 



misclassified close binaries (Miszalski 2009). Further exam- 
ples will be presented elsewhere and we expect further high 
resolution spectroscopic observations of so-called wels will 
identify further close binaries. 



3.2 Lightcurve 

Time series photometry of ETHOS 1 was obtained with 
the MEROPE camera (Davignon et al. 2004) on the 1.2- 
m semi-robotic Flemish Mercator Telescope (Raskin et al. 
2004) from 24 August 2009 to 4 September 2009. Figure 0] 
shows the lightcurve phased with our ephemeris: HJD (min. 
I) = 2455076.0312±0.0007 + (0.53512±0.00019) E. The IDS 
and FORS spectra taken earlier during the same year cov- 
ered phases 0.40-0.44 and 0.80-0.85, respectively. As for the 
Necklace, the observed periodic variability is clear evidence 
of binarity. A sine fit to the lightcurve gives a very high am- 
plitude of 0.816 mag that indicates a very hot primary irradi- 
ating a cool main-sequence companion (see e.g. Fig. 1 of De 
Marco et al. 2008) . Such a large amplitude places ETHOS 1 
amongst a small group of post-CE binaries showing extreme 
irradiation effects with only Sab 41 having a larger ampli- 
tude (I = 0.849 mag, Miszalski et al. 2009a). The next 
largest amplitudes are found in the Necklace (/ = 0.747 
mag, Corradi et al. 2010) and K 1-2 (V = 0.680 mag, Ex- 
ter al. 2003). There are no other post-CE binaries without 
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visible nebulae that come close to these extreme amplitudes 
simply because they have evolved significantly further along 
the white dwarf cooling track (Aungwerojwit et al. 2007; 
Shimansky et al. 2009). 



4 THE NEBULA 

4.1 Morphology and Kinematics 

Narrow-band images of 300 s duration each were obtained 
during the VLT visitor mode program 085.D-0629(A) on 
19 June 2010 with the FORS H_Alpha+83, OIII+50 and 
OII+44 filters. The central wavelengths and full widths 
at half maximum (FWHM) of each filter are 656.3/6.1 
nm, 500.1/5.7 nm and 377.6/6.5 nm, respectively. Note the 
H_Alpha+83 filter includes the [N II] A6548, A6583 emis- 
sion lines. Figure [5] presents a montage of the images which 
show three components to ETHOS 1: (i) The main nebula 
19.4" across measured at 10% of the peak Ha intensity, (ii) 
the jets measuring 62.6" tip-to-tip also at 10% of the peak 
Ha intensity, and (iii) a faint bipolar outflow best visible in 
[O III] (Fig. 

The absence of an inner ring of low-ionisation structures 
(LIS) sets ETHOS 1 apart from the other nebulae surround- 
ing high-amplitude close binaries including K 1-2 (Exter et 
al. 2003), Sab 41 (Miszalski et al. 2009b) and the Necklace 
(Corradi et al. 2010). As the UV radiation field is compara- 
bly high in the other objects this may indicate the ring never 
existed rather than recent ablation of an LIS ring. There are 
in fact no low-ionisation species detected in ETHOS 1 ex- 
cept for in the shocked jet tips described in Sec l4.2l A faint 
detection in the [O II] image is due to the nebula continuum 
since there is no [N II], [S II] or [O II] emission recorded by 
our INT spectrum of the main nebula (Fig. [2]). 

Spatially resolved long slit echelle spectra of ETHOS 1 
were made using the second Manchester Echelle Spectrom- 
eter on the 2.1-m San Pedro Martir Telescope (Meaburn et 
al. 2003). The primary spectral mode was used with nar- 
rowband filters to isolate the [O III] A5007, Ha and [N II] 
A6548, 6583 emission lines. Observations were performed on 
23 September 2010 with a 2' x 2" slit and the Thomson 
2Kx2K CCD which gave a spatial scale of 0.38" /pixel and a 
resolving power of R ~ 30000 (2.82 km/s per pixel). A series 
of 1800 s exposures were taken with the slit oriented approx- 
imately along the minor (PA=60°) and major (PA=147°) 
axes. 

Figure [6] presents the emission line position-velocity 
(PV) diagrams observed. Alongside these are the same PV 
diagrams extracted from a SHAPE (Steffen & Lopez 2006) 
best-fitting spatio-kinematic model of ETHOS 1 built fol- 
lowing Jones et al. (2010). Figure[7|shows the model that was 
inspired by the faint Hb 12-like bipolar extensions (Vaytet 
et al. 2009) that are visible in Fig. [B- The basic compo- 
nents are a spherical outflow and nested thin-walled bipolar 
outflows. The inner bipolar outflow helps reproduce the X- 
shape and the outer (possibly related) Hb 12-like bipolar 
outflow protrudes from the central sphere. Our approach 
was not tailored to each emission line observed, but instead 
aimed to create a good overall approximation to the features 
seen in Figs. [5] and [6] The jets were not incorporated in our 
model but they appear to be well-aligned to the major axis 




Figure 7. A SHAPE spatio-kinematic model of ETHOS 1. 

unlike other PNe (e.g. Hb4). The faint bipolar extensions 
(Fig. (B) were included in our model but their low surface 
brightness precluded their detection on the 2.1-m telescope 
at this resolution. We note that no attempt has been made 
to replicate the relative surface brightnesses of each compo- 
nent in the model and deeper observations should be able 
to detect these fainter features. 

Along the major axis the PV diagrams of the inner neb- 
ula display a figure-of-eight profile best seen in Ha that is 
typical of a bipolar outflow with a 'pinched' waist (e.g. Icke, 
Preston & Balick 1989). This bipolar profile is also appar- 
ent in [O III], but is far more irregular and lacks emission 
from the central region. The southern component splits into 
two velocity components that are consistent with the front 
and back walls of a hollow bipolar shell. A velocity ellipse is 
also seen in the [O III] profile surrounding the figure-of-eight 
profile. This is consistent with a spherical shell of material 
and the minor axis PV diagrams support this interpretation. 
Along the minor axis the outer ellipse represents the spheri- 
cal shell, while the inner knot at negative velocity represents 
the pinched bipolar shell (e.g. Vaytet et al. 2009). There is 
a small amount of asymmetry in the minor axis which is 
partly due to imperfect alignment of the slit along the axis 
and partly intrinsic. 

No emission is detected from the central region in [N II] 
consistent with our other observations. There is however an 
elongated excess of Ha emission in the inner nebula aligned 
along the major axis that is not matched in [O III] (e.g. Fig. 
(B)- While this could be explained by ionisation stratifica- 
tion (e.g. the same feature may be visible in O V or O VI), 
it may also be related to similar but 'tilted' inner nebulae in 
PN G126.6+01.3 (Mampaso et al. 2006), M 2-19 (Miszalski 
et al. 2009b) and A 41 (Miszalski et al. 2009b; Jones et al. 
2010). We suggest these nebulae could be the remnants of 
a precessing jet (Haro-Corzo et al. 2009; Raga et al. 2009) 
or nozzle (Balick 2000). The uncanny resemblance between 
PN G126.6+01.3 and the models of Haro-Corzo et al. (2009) 
supports this hypothesis. 

The symmetry axis of our best-fitting SHAPE model has 
an inclination of 60 ± 5 degrees to the line of sight. This 
inclination is consistent with the non-eclipsing lightcurve, 
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Table 2. Distance dependent kinematic ages of nebulae and jets 
of PNe with close binary nuclei. 



Name ^nebula %ets Reference 

(yrs/kpc) (yrs/kpc) 

A 63 3500 ± 200 5200 ± 1200 Mitchell et al. (2007) 

Necklace 1100 ± 100 2350 ± 450 Corradi ct al. (2010) 

ETHOS 1 900 ± 100 1750 ± 250 this work 



4.2 Emission line diagnostics 

Table [3] gives our measured emission line intensities for our 
extracted IDS and FORS spectra. Dereddened intensities 
were calculated using the Howarth (1983) reddening law and 
c = 0.18. As Exter et al. (2003) point out an abundance 
analysis for such high excitation nebulae requires photoion- 
isation modelling to properly measure the chemical abun- 
dances. Since our spectral information is rather limited and 
lacks reliable UV and far-blue coverage we postpone this 
work for now. Abundances of the jets are likely heavily dis- 
torted by their shocked nature, but we can investigate their 
emission line ratios for diagnostic purposes. 

The nebular package in iraf (Shaw & Dufour 1995) 
was used to determine electron temperatures and densities 
in the main nebula and SE jet. Assuming n e = 1000 cm -3 , 
the main nebula has a very hot T e = 17700 ± 500 K. This is 
comparable to K 1-2 with T e = 17500 ± 1000 K, but is larger 
than T e = 14800^4^ K measured in the Necklace (Corradi 
et al. 2010). The [Ar IV] A4711/4740 density diagnostic from 
our VLT spectrum lies near the lower density limit so we find 
n e — 850 ± 1000 cm -3 . We also find the jets to be cooler 
T e = 12900 ± 1000 K and to have a similarly low density by 
using the [S II] A6731/6717 diagnostic to provide an upper 
limit of n e < 1000 cm" 3 . 

The emission line ratios from separate extractions of 
the jet tips in Tab. [3] can be compared to model predic- 
tions (Dopita 1997; Raga et al. 2008). These works aim to 
reproduce FLIERS as the result of shocks in the strongly 
photoionised medium provided by the CSPN giving the ap- 
perance of a 'reverse-shock'. In ETHOS 1 the presence of 
[O III] in the jet and its relatively low T e =12900 K indi- 
cates strong radiative cooling is occurring in the post-shock 



the high degree of mirror symmetry evident in the main 
nebula about the minor axis and the large visibility fraction 
of the jets. A heliocentric radial velocity of — 20±5 km/s was 
determined from our model for the nebula. The expansion 
of the model nebula follows a Hubble-type flow, where all 
velocities are proportional to the radial distance from the 
nebular centre, equivalent to an expansion velocity of 55 
km/s at the edge of the spherical component. This is above 
average for a PN but not unexpected for a bipolar outflow 
(Corradi & Schwarz 1995) with He II A6560 present (Richer 
et al. 2008). Figure [3] shows He II A4859 and we confirm the 
detection of He II A6560 in our MES spectroscopy. 

Although the distance to ETHOS 1 is unknown, the 
spatio-kinematic model was used to determine the kinemat- 
ical age per kpc of the inner nebula as 900±100 years kpc -1 . 
Similarly, assuming the jets have the same inclination as the 
internal bipolar structure, we find their kinematical age to 
be 1750±250 years kpc -1 at their tips, almost double that of 
the central region! The velocities of the jet tips reach — 55±5 
km/s (SE jet) and 65 ±5 km/s (NW jet). Their deprojected 
velocities at our assumed inclination are —110 ± 10 km/s 
and 130 ± 10 km/s, respectively. 
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Figure 5. A montage of VLT FORS ETHOS 1 images, (a) Colour-composite of Ha+[N II] (red), [O III] (green) and [O II] (blue); (b) 
and (e) Ho+[N II]; (c) and (f) [O III]; (d) [O II]; Each image measures lxl arcmin 2 with North up and East to the left. 
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Figure 6. Model and observed position-velocity diagrams for the major (PA=147°) and minor (PA=60°) axes of ETHOS 1. The velocity 
scale is centred on the model-derived heliocentric radial velocity of —20 km/s. The jets are not incorporated into the models and our 
observations are not sensitive to the faint bipolar outflow (major axis model). The orientation of the major and minor axes are North is 
up and East is up, respectively. 
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Table 3. Measured (-Fa) and dereddened nebula emission line intensities using c = 0.18. 



Component: 


Main Nebula 


SE Jet Tip 


NW Jet Tip 


Main Nebula 


SE Jet Whole 


Instrument: 




IDS 




IDS 




IDS 


FORS 


FORS 


Identification 




h 






Fx 


h 


Fx 




F A 




[in e v j o*±zo 1 


122 


140 


- 


- 


- 


- 


- 


- 


- 


- 


[W 11J O ( £i I ' 


- 


- 


642 


713 


133 


148 


- 


- 


- 


- 


\Nr TTTl ^SfiQt 


38 


42 


133 


147 


67 


73 


- 


- 


- 


- 


TNTrv TTTl J-T-T^T /TT QQ^ftt 
[INC 111J -)-rlCl/ 11 oyoo ' 


38 


41 


- 


- 


- 


- 


- 


- 


- 


- 


T-TX A1 01 t 

no 4iui 1 


40 


44 


- 


- 


- 


- 


25.8 


27.8 


28.8 


31.0 


nc ii 4zuu 


- 


- 


- 


- 


- 


- 


4.3 


4.6 


- 


- 


TT AOdT 


- 


- 


- 


- 


- 


- 


3.2 


3. 1 


- 


- 




51 


54 


- 


- 


- 


- 


45.7 


48.1 


45.7 


48.2 


U ±± J-J lOUO 


* 


* 


- 


- 


- 


- 


9.7 


10.3 


11.1: 


11.7: 


T-To T AA71 
ne 1 / 1 


- 


- 


- 


- 


- 


- 


- 


- 


6.4: 


6.6: 


nc ii QuQi 


- 


- 


- 


- 


- 


- 


3.9 


4.0 


- 


- 


nc ii 40ou 


115 


117 


20: 


21: 


50: 


51: 


111.8 


113.8 


61.1 


62.2 


TAr TVl 471 1 
[-rt-I 1 V J t i 1 1 


* 


* 


- 


- 


- 


- 


11.8 


12.0 


- 


- 


PWo TVl 479^ 
[IN C 1 V J ^± ( ZO 


- 


- 


- 


- 


- 


- 


2.2 


2.2 


- 


- 


f A -r TVl A7AC\ 


* 


* 


- 


- 


- 


- 


9.0 


9.1 






T-T R A Qdl 

np 4ooi 


100 


100 


100 


100 


100 


100 


100.0 


100.0 


100.0 


100.0 


TO TTTl AQ^Q 
iu iiij ^y<jy 


128 


127 


317 


313 


317 


313 


132.5 


131.2 


299.5 


296.6 


[O III] 5007 


404 


398 


975 


960 


960 


946 


393.1 


387.3 


891.5 


878.3 


He II 5412 


10: 


10: 




_ 


_ 


_ 


8.9 


8.4 


4.2: 


4.0: 


[O I] 6300 


_ 


_ 


■ 


* 


33: 


30: 


- 


- 


- 


- 


[Ar V] 6434 


5: 


1: 


_ 


_ 


_ 


_ 


_ 


_ 


_ 


_ 


[N II] 6548 


_ 


_ 


57 


50 


79 


69 


- 


- 


- 


- 


Ha 6563 


317 


278 


304 


266 


293 


257 


_ 


_ 


_ 


_ 


[N II] 6583 


_ 


_ 


204 


179 


273 


239 


- 


- 


- 


- 


[S II] 6717 


- 


- 


* 


* 


18: 


12: 










[S II] 6731 






* 


* 


43: 


38: 










[Ar VI 7005 


12: 


10: 


















[Ar III] 7135 


11: 


9: 


37 


31 


12 


36 










[Ar IVJ+He II 7175 


* 


* 


















N I 7452 






59 


■19 














log([0 III] 5007/Ha) 




0.16 




0.56 




0.57 










logQO III] 5007/H/3) 




0.60 




0.98 




0.98 










logQN II] 6583/Ha) 








-0.17 




-0.03 










log([S 11} /Ha) 












-0.50 










log([S II] 6731/Ha) 












-0.83 










logQO I] 6300/Ha) 












-0.94 










logQO II] 3727/ [O III] 5007) 








-0.13 




-0.80 










logQO II] 3727/H/3) 








0.85 




0.17 











^Lines more uncertain in IDS spectra because of bright sky background. 
:Low S/N measurement. 

*Line present but too low S/N to be measured. 



region. We find the line ratios agree well with log Ui ~ —4.4 
from Dopita (1997) and the measured values of FLIERS in 
NGC 6543 and NGC 7009 (Raga et al. 2008). The Raga et 
al. (2008) models systematically underestimate the [O III] 
A5007/Ho: ratio, although this is probably due to the rel- 
atively low stellar temperatures adopted in their models. 
The spatial distribution of the key emission lines Ha, [N II] 
A6583 and [O III] A5007 are also consistent with the Raga 
et al. (2008) models. Figure [8] presents spatial profiles of 
these lines where the CSPN continuum was subtracted from 
extractions adjacent to the lines. The results qualitatively 
agree with both IC 4634 and the models presented by Raga 
et al. (2008) with [N II] extending ~2" further out than the 
other emission lines. This effect is also evident in Fig. 



5 CONCLUSIONS 

We have introduced ETHOS 1 (PN G068. 1+11.1), the first 
spectroscopically confirmed PN discovered from the Ex- 
tremely Turquoise Halo Object Survey (ETHOS, Miszalski 
et al. in prep). An irradiated close binary central star was 
discovered with an orbital period of 0.535 days and an am- 
plitude of 0.816 mag. The extreme amplitude is second only 
to Sab 41 (Miszalski et al. 2009a) and is consistent with the 
presence of a very hot central star that produces the highly 
ionised nebula (T e = 17700 K). The Necklace (Corradi et 
al. 2010) and K 1-2 (Exter et al. 2003) also share similarly 
large amplitudes and high-ionisation nebulae, although the 
absence of low-ionisation filaments in ETHOS 1 may suggest 
a slightly different evolutionary history. VLT FORS spec- 
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Spatial axis (arcsec) 

Figure 8. Spatial profile of the Ha, [N II] A6583 and [O III] 
A5007 emission lines from the INT IDS spectrum. The insets show 
normalised profiles of the ends of the jets. No dereddening has 
been applied. 



troscopy of the CSPN confirms the presence of a close bi- 
nary with a large velocity separation between primary and 
secondary components and the nebula. The presence of N 
III, C III and C IV emission lines continues the trend seen 
in other irradiated close binary CSPN. These weak emission 
lines are typical of many CSPN classified as wels in the liter- 
ature and we expect many of these will turn out to be close 
binaries. Further observations are required to constrain the 
orbital inclination, masses and radii of the binary CSPN. 

A spectacular pair of jets travelling at 120 ± 10 km/s 
accompanies the inner nebula of ETHOS 1 adding further 
evidence towards the long suspected relationship between 
binary CSPN and jets. Their tips present emission line ra- 
tios that are consistent with shocked models of FLIERS in 
PNe. The fact that the jets are detached suggests a limited 
period of jet activity consistent with a transient accretion 
disk before the CE is ejected. The kinematic age of the jets 
(1750 ± 250 yrs/kpc) was found to be older than the in- 
ner nebula (900 ± 100 yrs/kpc) supporting this hypothesis. 
ETHOS 1 therefore continues to follow this trend previously 
identified in A 63 (Mitchell et al. 2007) and the Necklace 
(Corradi et al. 2010). Both ETHOS 1 and the Necklace have 
younger kinematic ages than A 63 consistent with the appar- 
ent ongoing cooling in the younger jets via [O III] emission. 
A close binary engine powering jets is out of place in the 
literature with most models of jets incorporating orbital pe- 
riods of several years (Cliffe et al. 1995; Haro-Corzo et al. 
2009; Raga et al. 2009). These models may be adequate if the 
accretion disk is established before the companion begins its 
in-spiral phase as our kinematic ages suggest. Nevertheless, 
new models with shorter orbital periods may be necessary to 
model the jets in objects like ETHOS 1 which are relatively 
more collimated than extant model jets. 

VLT FORS imaging and MES high-resolution spec- 
troscopy of the inner nebula of ETHOS 1 was conducted. 
The A-shaped inner nebula was reconstructed using a SHAPE 
kinematic model consisting of bipolar outflows and a spher- 
ical nebula. The faint bipolar extensions particularly visible 
in the [O III] FORS image could be an earlier ejection oc- 
curring at about the same time as the jets or alternatively 



may represent a bubble produced by the jets clearing out a 
cavity around the PN. Similar bipolar outflows are seen in 
A 65 (Walsh & Walton 1996), PPA 1759-2834 (Miszalski et 
al. 2009b) and Fg 1 (Boffin & Miszalski in prep.) suggesting 
they may have a binary-related origin. It may also be possi- 
ble that ETHOS 1 represents a slightly more evolved state 
of Hb 12 (Vaytet et al. 2009) where the spherical nebula is 
ejected on top of a pre-existing Hb 12 system. 
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